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Abstract 

We derive a relation between leading finite size corrections for a 1+1 dimensional 
quantum field theory on a strip and scattering data, which is very similar in spirit 
to the approach pioneered by Liischer for periodic boundary conditions. The consis- 
tency of the results is tested both analytically and numerically using thermodynamic 
Bethe Ansatz, Destri-de Vega nonlinear integral equation and classical field theory 
techniques. We present strong evidence that the relation between the boundary state 
and the reflection factor one-particle couplings, noticed earlier by Dorey et al. in the 
case of the Lee- Yang model extends to any boundary quantum field theory in 1+1 
dimensions. 



1 Introduction 

Finite size corrections are important tools in QFT and their use is mandatory in interpret- 
ing numerical data. Nearly two decades ago M. Liischer [T] gave a general expression for 
the leading finite size correction of particle masses and scattering states valid in any QFT 
with periodic boundary conditions, in terms of the infinite volume scattering data. 
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Liischer's result was specialized to 1+1 dimensional field theory 0, using the exact 
scattering amplitudes known in integrable field theories. For example, in a theory with a 
single particle of mass m and elastic scattering amplitude S (6) the leading correction to 
the particle mass in finite volume is the following: 




(similar formulae exist for general particle spectra). The first term (which, if present, 
dominates over the second) appears whenever the particle occurs as a bound state of itself, 
and is proportional to the residue of the S'-matrix at the pole, while the second term is 
present in general and is of order L~^/^e~'"^. 

In integrable 1+1 dimensional models restricted to a spacetime strip of finite width L 
with integrable boundary conditions, some effort has been devoted recently to the deter- 
mination of the leading corrections to low lying - mainly ground-state - energy levels as 
functions of L in classical field theory [Sill]. However, these investigations did not look 
for a general, model independent form of the asymptotic correction to the infinite volume 
ground state energy, which is the goal of the present work to consider. In addition, the 
results of the present work will be derived at the quantum level, and it is shown that they 
agree to the results in [Hllll in the classical limit. 

In this paper we derive a formula, similar in spirit to Liischer's result in the case of 
periodic boundary conditions, which expresses the finite size correction to the energy of 
a state containing no particles (but possibly excited boundary states) in terms of the re- 
fiection amplitudes on the semi-infinite line. Once again, the leading term depends on 
the presence of certain poles in the refiection factors. When such poles are present we 
get corrections of the form e"*"^ (for a particle of mass m) with coefficients that depend 
on the residue of the respective pole, otherwise the leading correction is given by a more 
complicated integral containing the refiection amplitudes and is of order L~^/^e~^™^. We 
show that these results are consistent with the boundary version of thermodynamic Bethe 
Ansatz (TBA). Using the general results we determine the form of leading finite size cor- 
rections in boundary sine-Gordon theory, which we check against results from a nonlinear 
integral equation (NLIE) in the Dirichlet case, and compare them to the classical field 
theory results in 011]. 

2 Cluster expansion and finite size corrections 

For simplicity let us consider a theory with only one scalar particle in the bulk. Assume 
further that the theory is defined on a strip of width L, with boundary conditions labeled 
by a and /? such that the elastic refiection amplitudes of the particle on these boundaries 
on the semi-infinite line, Ra (9) and Rp (9) are known. To derive the asymptotic (L oo) 
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form of the ground state energy £'°^(L), we consider the partition function of the theory 
on a cylinder in two alternative forms: 

Z^p {R, L) = Tre-^^'^''^^) = («|e-^^(^) 

where H{R) denotes the Hamiltonian of the theory in the crossed channel (on a circle 
of circumference R), \a) and are the appropriate boundary states a la Ghoshal and 
Zamolodchikov 0, while Hap{L) is the Hamiltonian on the strip. In the limit i? ^ oo we 
obtain 

The right hand side can be expanded for large L in a so-called cluster expansion (the 
analogue of which was recently used in [^1 for the TEA in the other channel to settle long 
standing issues concerning the boundary entropy function). Inserting a complete set of 
eigenstates of the Hamiltonian H{R) on the circle we obtain 



(a 


\n){n\[3) 


(n 


\n) 



"S'^^bulk 

For our purposes, only the first few terms are interesting: 

^ ' (0|0) 

{a\9n){0n\(^) ^-L{Eo {R)+m cosh 6»„) 



{On\On) 

n 



-L{Eo{R)+m cosh 6„+m cosh 9^) 



+ three particle contributions + .. . (2-1) 
In the limit R ^ oo the boundary state |/?) has the general form 

\f3) = (l + (^) (-^) (^) + • • •) |0) (2-2) 

(similarly for where the dots denote terms containing more particles, (9) creates an 
asymptotic particle of rapidity 9 and 

K,{9) = Rfs(^-^-9 

Normalizing the bulk vacuum energy to and collecting the leading (two-particle) terms 
for large L yields 

poo jn 

Z^p{R,L) = 1 + mR —cosh9Ka{9)Kp{9)e~^"'^'°'^^ + ... 

Jo 27r 



K^i9) = K^(-9) = R^[- + 9 



where we replaced the sums by integrals over the density of states 

> niR I — cosh 

^ J 271 

On 

and used real analyticity of the reflection factor 

R{e)* = R{-e*) 

For the energy we obtain 

/■oo jn 

E^p{L) = -m — cosh^ir,(^)K^(^)e-2™^=-^^ (2.3) 

One can generalize the above argument to the case when there is a multiplet of particle 
of mass m. In that case the reflection factor is a matrix 

id): 

where r and s denote the multiplet labels. The boundary state takes the form 



/•CO jn \ 



where 

(ey = - e 

and s denotes the charge conjugate of multiplet member s. The leading finite size correction 
takes the form 

E^p{L) = -m ^ cosh^Tr {K^ (9) {9)) g-^-^-^^^ (2.4) 

Let us turn now again to the case of diagonal scattering. If the theory is integrable 
then the boundary state is known exactly 

1/3) = exp (^j^ K, (9) At i-9) (9) ^ |0) (2.5) 

and summing up the cluster expansion one can derive a TBA equation for the ground state 
energy as done by Leclair et al. in 0, which is analyzed in detail in the next section. The 
large volume limit of this TBA is exactly given by the formula (j2.3|) . 

However, if there is a pole in the refiection factor at the imaginary rapidity 9 = in/ 2 
of the form 

R^O) ^ ^ (2.6) 
the boundary state has an extra one-particle contribution 
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/oo JO \ 
-K,{e)AH-e)AHe) + ...j\o) (2.7) 

or in the case of an integrable field theory one has 

(poo JQ \ 
X ^^/3WAtM)At(^^)j|0) (2.8) 

(due to the exclusion principle, adding the one-particle term to the exponent makes no 
difference). Ghoshal and Zamolodchikov [HI identify gp with gp, but Dorey et al. in |H1 
found the relation 

9. = f (2.9) 

comparing one-point functions calculated from the boundary state using form factor meth- 
ods to TCSA results in the case of the Lee- Yang model. As we show below, our results are 
fully consistent with this latter suggestion. Namely, using the state ()2.7|) in the expansion 
p. HI one gets for the leading term the result 

Z^p {R, L) = l+ mR-g^-gpe-^'^ + ... 

and so 

E^p (L) = -m^e-"^^ + . . . (2.10) 

The existence of this contribution requires that both left and right boundary reflection 
matrices have a pole at ^ = 111/2. It is also important to note that if one-particle terms are 
present in the boundary state, the two-particle term of the cluster expansion is divergent 
and needs to be regularized. 

Hereafter we shall refer to ()2.3|2.in|l as the boundary Luscher formulae because they 
are a natural generalization of Liischer's original results ^ to the case of quantum field 
theory defined on a strip. 



3 BTBA in the infrared limit 

In order to lend support to relation ()2.9p between the one-particle coupling constant ap- 
pearing in the refiection factor and the coefficient of the one-particle term in the boundary 
state, we calculate the infrared limit of boundary TBA (BTBA). Although the original 
derivation given by Leclair et al. [Z| only considers the case when there is no one-particle 
term in the boundary state, it can be easily argued that the presence of one-particle terms 
makes no difference to the end result, and this is also supported by numerical studies using 
comparison with truncated conformal space (TCS) by Dorey et al. in pj. 

For simplicity we consider a theory with a single particle of mass m on a strip of length 
L. The BTBA equation takes the form 

poo jnl 

e{e) =2mLcoshe - {6 - e')\og (l + x{0') e-^'^'n (3.1) 
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where 

is the derivative of the phase of the two-particle S'-matrix S{9), and 

using the notations of the previous section. Using crossing-unitarity (H] 

K^{9) = Si29)K^{0) (3.2) 

and unitarity 

S{9)Si-9) = 1 
it is easy to see that x(^) is an even function. 



The energy of the ground state is given by 

/oo jn 
— cosh e log {l + x{9)e'<'^) 



(3.3) 



In the case where there is no one-particle coupling, x(^) is regular on the real axis and for 
large L 

e (9) = 2mL cosh 9 + (e'^"*^) (3.4) 
which gives the following asymptotics for the energy: 

/oo jn 
— cosh {9) ^ee-^""^'"'^^ + O (e-^™^) 
■oo 47r 

in perfect agreement with ((231) • The same is true for the case when only one of the 
boundaries has a nonzero one-particle coupling. The corresponding K factor has a simple 
pole at ^ = 0, but since the other K is regular and x is an even function, x has no 
singularity at all. This can be understood by recalling that generally 

S{0) = -1 

so crossing-unitarity (13. 2|] entails 

KaiO) = -K^{0) 

and as a result, K factors must have either a zero or a pole at 6^ = 0. Even if there is 
a first order pole in one of the K factors corresponding to a one-particle coupling to the 
boundary, the product of the two K factors is still regular. 

However, for a theory with nonzero one-particle coupling on both boundaries x has a 
second-order pole at 6* = 0. The logarithmic terms in ()3.H) and ()3.3|) remain integrable and 
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the BTBA makes perfectly good sense, but to obtain the correct asymptotic expansion one 
needs to be careful. For large L, the energy can be written in the form 



/OO J^Q 
— cosh ^ log + e" 



2m,L cosh ( 



) 



(3.5) 



Now we cannot use log(l + x) ~ x to expand the logarithm as x(6') diverges at 6* = 0. In 
fact, using ()2.6|1 this divergence can be computed as 



X{0) 



dig} aWp 



402 4sinh2 
We can use the following integral formula [H] 



dx log 



to evaluate the energy in the following way 



2Ti{a-h) ,a,b>0 



(3.6) 



Ea/3{L) 



-m I — cosh tl 



An 



log 1 + 



gig} 
^( 

4 sinh 9 



-2mL 



+ log 



X{d)e 



'2mL cosh 6 



1 I -2mL 
^ 4sinh^e^ 



igagpl -rnL 

-m — ■ — e — m 



— cosh 9 



log 



1 + X{0)e 



-2mL cosh 9 



1 + 



-2mL 



(3.7) 



The remaining integral is an expression of order e 2*"^, while the leading term agrees with 

^MM for gagp > o. 

It was already noted in (H] (using a comparison with truncated conformal space method) 
that the BTBA equation only gives the correct ground state energy in this case, and for 
values of boundary parameters such that gag/s < a suitable analytic continuation must be 
applied.^ Note also that in the definition of the boundary coupling from ()2.6j) a branch 
of the square root function must be chosen. In all known cases (e.g. Lee- Yang in |9| or 
sine-Gordon in the following section) there exists a branch choice such that the boundary 
couplings depend analytically on the boundary parameters. In this case a straightforward 
analytic continuation from the range of parameters for which gag/s > gives the result 

^ai3{J^) — —rn——e +... 



in the parameter range where gagp < 0, which is consistent with ()2.9I2.1()|) . Furthermore, 
the result ()3.7|] also yields a formula for the regularized two-particle contribution. The two- 
particle contribution is regularized by the logarithm in ((^1^1), which represents a partial 

^ Using the notation of paper 0, the regime gagp > corresponds to the boundary parameter range 
—2<b<2 (for boundary conditions a = I and /3 = ^{h), with the boundary parameter h related to 6), 
which is the range in which they found that BTBA compares well with truncated conformal space. For 
parameters out of this range the authors needed to resort to analytic continuation in order to match the 
numerical data. 
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re-summation of the higher multi-particle contributions, and the bulk interaction plays 
no role at leading order in large volume. We conjecture that the leading asymptotics 
calculated from (|3.7I1 is valid also for non-integrable theories, since for large L the leading 
contribution comes from a very small region around 6 = Q and for such small values of the 
rapidity the inelastic scattering channels should play no role provided m is the smallest 
particle mass in the spectrum. 



4 Application to sine-Gordon theory with integrable bound- 
ary conditions 

Next we apply these results to sine-Gordon theory on a strip with integrable boundary 
conditions of Ghoshal-Zamolodchikov type IS]- The model has a complicated spectrum 
with non-diagonal bulk and boundary scattering. However, to derive the leading term 
p.3|2.in|l we used only ()2.2I2.7|1 which is valid in general field theories (even non-integrable 
ones) in two dimensions, as follows from the derivation of (H]. Integrability is only needed 
to derive the exact exponentiated form of the boundary state (j2.5l2.8jl . Besides that, the 
derivation of Section 2 carries over to the case of non-diagonal boundary scattering in a 
straightforward manner, while the effects of bulk scattering appear only in higher correction 
terms. 

The action of the model is 



A 



/OO r nh 

dt / dx/:5G(x,t)-l4°^($(0,t))-Vi^)($(L,t))-a 



d<^{L,t) 
dt 



(4.1) 



where the bulk Lagrangian density is given by 



Csg{x, t) = - (9^$(x, t)Y- — {l- cos (3^) 
and the boundary interaction is described by the potentials 



l-cos^($(2;,t) 



0, L 



(4.2) 



(the Mz oo limit corresponds to Dirichlet boundary conditions). The last term (pro- 
portional to a) is allowed by integrability, and on the strip - in contrast to the theory on 
the half line - it cannot be eliminated from the action by an appropriate "gauge" transfor- 
mation [in]. The spectrum is constituted by the solitons and breathers, and the scattering 
theory on the half line is known. The reflection amplitudes on the two boundaries can be 
parameterized with the Ghoshal-Zamolodchikov parameters rj^ and {z = 0, L). 

Using a perturbed conformal field theory description one can write the action of the 
quantum theory as 



A 



pCFT 



/oo /"O /"OO n 

dt j dx : cos /?$ : +/io / dt : cos - ($(0, t) - $c 
■oo J —oo J —oo 2 
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oo 



+ fiL I rft: cos^($(L,t) : 



J — oo J —oo ^ J —oo 



dt 



The parameters of this action are related to the boundary parameters ?7, d by the UV-IR 
relation [TT] 



sm 



(-1)"/^ sin— ^sinh- 



/icrit 2 A + 1 A + 1 

P^Z P^Z Vz . I^Z /X 

COS = cos cosh (4.3) 



/icrit 2 A + 1 A + 1 



where 



flcrit = J^^ , A = ^-l (4.4) 

sm ^ P 



and the sign factor —( — 1)^/^ arises from the fact that on the left end {z = 0) of the 
strip the role of the soliton and antisoliton are interchanged with respect to the right end 
{z = L). In the classical limit flz can be replaced by and fi by w? / [3'^ from (I4.1I1 . and 
the 7], d parameters must be scaled as JHl 

77 = r/,i(l + A) , ^9 = ^ci(l + A) (4.5) 

keeping r/cb ^ci finite. This leaves us with the classical UV-IR relations 

■ t ^^z/L . ■ . a 

sm =— (—1)' sm?7cismhWci 



Merit 2 



^ cos ^ = cos T^ci cosh T^cl (4.6) 



M„it 2 

4m 



Merit 



In the regime < 47r the spectrum contains the first breather, which has the refiection 
matrix [12] 



'^^^^ (^ + 1) + + ' sinh(l-.f) 



The breather refiection amplitudes are all independent of the parameter a, which appears 
only in the soliton refiection matrices. The refiection factor ()4.7|l does have a pole at 
9 = m /2 with the singular term taking the form 

r^a),n^ . .1 + COS - sin ^ „ ?7 2 -(9 1 

B}nHe) ~ Ai- tan^ -V tanh^ — — 

\^ ^ ' 1 - cos ^ + s n 1^ 2A 2A 29 - tir 



which gives the following result for the coupling of the first breather to the boundary with 
parameters rj and 'd 



/ 1 + cos 1^ — sin 7^ n ^ -d , , 

gi(r],m=2j ^ ^ tan -V tanh— 4.8 

^ y 1 - cos § + sin 1^ 2A 2A ^ ' 

(The expression under the square root is always positive as long as A > 1 which is necessary 
for the first breather to exist in the spectrum). However, some care must be taken, because 
the sign of the coupling gi must be opposite on the two ends of the strip since these are 
related by a spatial reflection under which the first breather is odd. 

As a result, formula ()2.inp together with ()2.9|) predicts that in the regime < in 
(A > 1) for generic boundary conditions the leading finite size correction to the ground 
state on the strip is given by 

EaAL) = -^mi(?i(r/z.,^L)(-^7i(r/o,^o))e-™^^ (4.9) 

1 + cos - sin 770 , ''^o Vl ^ m t 

= mi ^ ^ tan tanh tan tanh ^e'^^^ + . . . 

l-cos§ + sin|^ 2A 2A 2A 2A 

where ?77i is the mass of the first breather, r^OjT^o and rji^'&L are the boundary parameters 
characterizing the boundary at the left and right ends of the strip, respectively. Soliton 
corrections are always subleading, since there is no one-particle coupling of solitons to the 
boundary and rui < 2M. 

It is also interesting to compute the one-particle term contributed by the second 
breather, since in the regime where the second breather exists (/5^ < Stt/S i.e. A > 2), 
its mass m2 = 2mi cos ^ is lower than 2mi and so it dominates over the two-particle 
correction from the first breather, giving the next-to-leading finite size correction. Fur- 
thermore, if any of the rj^ (-2 = 0, L) parameters is zero, then the leading term ()4.9|1 
vanishes, leaving the second breather's one-particle term as the leading one. 

The refiection amplitude for the second breather is of the form [T2j 



(I) (^ + 1) + il.) 



(v^ _ I L\ (ii^ _ 1 _ J-^i (IK - I ^ l.\ (il^ _ 1 _L J."! 

VttA 2 2A/ V ttA 2 2A/ WA 2 2A/ V ttA 2 2A/ 

(IK - l.\ (i^ 4-1 - ±\ (IK A_ 1 A_ ±\ (r&z I i I J_^ 

VttA 2 2AM ttA 2 2A/ KttX 2 2A M ttA 2 2A / 

from which the boundary coupling can be extracted 



(4.10) 



^g) _ 2 tan (jL - ^) tan + ^) tan [f^ - g) tan [f^ + g) 



tan^^tan|^tan(f + ^) 
and the one-particle contribution to the finite size corrections becomes 

_ ^^^2(^0,^0) ^2 (r^L,^.)^_^,, (4.12) 
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(here the sign of the couphngs are the same at the two boundaries as the second breather is 
even). We remark that in this model there is in fact no sign ambiguity in this couphng once 
the signs of bulk three-particle couplings are fixed, because the soliton reflection factors 
contain a pole at ^ = n = 1,2, . . . ,[X\ whose residue is proportional to 

where is the bulk soliton-antisoliton-breather coupling, while gn ij], -d) is the one- 
particle coupling of the nth breather to the boundary. As gn ij], -d) appears linearly, it 
can be unambiguously calculated and the result matches the formulae (14.814.111) . We note 
that in the limit 'd ^ oo (Dirichlet boundary conditions) the couplings simplify to 



'l + cos§ -sin|^^^_^ ^ 



giiri, oo) = 2 J ^ ^tan- 

l-cos^ + sin|^ 2A 

2tan(^-:^)tan(^ + ^) 
g2{r],oo) = (4.13) 



tan^L^tan^tand + 1^) 

5 Comparison with the Dirichlet NLIE 

In this section we check the boundary Liischer formulae derived previously by comparing 
them to the exact ground state energies at large but flnite L-s in the sine-Gordon model 
with Dirichlet boundary conditions. We determine the exact ground state energies from 
the (ground state version of the) NLIE proposed recently in [15], which generalizes [Z| by 
allowing the sine-Gordon field to take two different values at the two boundaries of the 
strip. 

The NLIE is a complex nonlinear integral equation of the form 

/oo 
dxG{e -x-i7]) log (l - e'^(^+^'')) 
-oo 

/oo 
dxG{e -x + ir]) log (1 - e^^^^-^")) (5.1) 
■oo 



where 



P{9\Ho,Hl) = 271 dx{F{x,Ho) + F{x,HL) + G{x) + J{x)) 



G{e) = / — e 



271 sinh ^ k cosh | k 



m = / T^e 



sinh ( ^^^TT^A; I cosh ( ^^^tt^A; 



27r sinh j^k cosh |/c 



^ 27T 2 sinh j^k cosh |/c 
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M is the soliton mass, and A is the parameter defined in ()4.4|1 . The parameters i^o.L can 
be expressed from the boundary values of the sine-Gordon field $: 

Ho = , H, = 

and the equation in the present form is valid only for < ifo,L < 1 + A^^. 

This equation must be solved for the unknown function Z{6). All the functions in the 
equations are analytic in some strip containing the real 9 axis, and r] must be chosen to lie 
inside that strip. Once Z{9) is obtained (usually by a straightforward iteration procedure), 
the ground state energy of the theory can be calculated using 

/°° dr 
— sinh(x + iri) log (l - e^^(^+^'')) (5.2) 
oo 27r 

It satisfies 

E{L) as L ^ oo 

which means that the vacuum energy extracted from the NLIE is normalized to zero 
vacuum energy density and zero boundary energy in infinite volume. The true vacuum 
energy (calculated e.g. in the perturbed conformal field theory framework) has the form 



EpcT — ^hulkL + ta + e/3 + E] 



NLIE 



where the bulk energy constant £^buik and the boundary energy terms ea,i3 are known exactly 
[TT] . The NLIE energy ()5.2|] is normalized the same way as the boundary Liischer formulae, 
so they can be compared directly. 

In the repulsive regime A < 1 one can calculate the leading large volume asymptotics 
by substituting 

Z {9) = 2ML sinh 9 + P{9\Ho, Hl) 
into ()5.2p and taking r/ ^ |. The result is 

/oo jn 
— cosh 9 (irf {-9) Kf {9) + Kl' [-9) Kf {9)) e-2*^^™«he ^5 3^ 
■00 

where 

Kf{9) = R^{i'--9y^ , Kf{9) = Rp{i'--9y_^ 

in terms of the soliton and antisoliton refiection factors found by Ghoshal and Zamolod- 
chikov 0. We made use of the relation [7| fTK] 

P{9\Ho,HL)=t log 5++ [29) - I log [9)1 -i log Rp {9)1 + n 

(where 5++ is the soliton-soliton phase shift) and also of bulk crossing symmetry and 
boundary crossing-unitarity to bring the result to the form (15. This is in perfect agree- 
ment with the result ()2.4|1 obtained from the cluster expansion. 
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Figure 5.1: Comparison of the Dirichlet NLIE with the boundary one-particle contributions 

The above procedure does not work in the attractive regime because the analyticity of 
Z does not extend far enough to allow taking the limit 77 — > |. In principle one should be 
able to recover the cluster expansion analytically from the NLIE expression (15 .21) even in 
the attractive regime; however, this requires more detailed understanding of the analytic 
structure of the NLIE ()5.H1 than what we have at present, and is out of the scope of the 
present paper. 

Instead we performed numerical comparison of NLIE results to the boundary Liischer 
formulae for several values of the coupling and the boundary parameters. An example is 
shown in figure l5?I] below, where A = 4 and so both the first and the second breathers are 
present in the bulk spectrum. Varying the boundary values of the sine-Gordon field we 
investigated the cases when one expects the asymptotics be given by the first or second 
breathers separately. The results are summarized in figureEHl where log(|£^|/M) is plotted 
against the dimensionless / = ML with M being the soliton mass. (The sign of the exact 
ground state energy always matched the sign of the appropriate asymptotic correction as 
calculated from ((T!9|) and ()4.12|l V 

For simplicity, we parameterize the boundary values $2, -2 = 0, L at the two ends of 
the strip in the following way 




and so ri^ = — (— l)^/^7r(l + A)</)2. The points marked by empty boxes, +-s, x-s and dotted 
circles are the exact ground state values computed from the NLIE when the boundary values 
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of the sine-Gordon fields on the upper and lower ends of the strip are 0o = — 0.170^ = 0.1, 
00 = 0.017 0L = 0.1, 00 = 0L = 0, and 0o = —0.08 0l = 0, respectively. The two 
upper continuous lines are the predictions of the first breather's asymptotic contribution 
(I4.9|l (with -(9 — > cxD to describe the Dirichlet case), while the two lower continuous lines are 
the predictions coming from the second breather's asymptotic contribution ()4.12|) when at 
one end of the strip rji = 0: 

The agreement with the NLIE data strongly supports the formulae derived from the clus- 
ter expansion. Although for the graphical comparison it is enough to plot only the second 
breather one-particle term, in order to make the agreement with the NLIE numerically 
satisfactory the two-particle contribution from the first breather also has to be included. 
The detailed numerical comparison is performed in the appendix and an excellent agree- 
ment is found. There's no such problem with the first breather term, however, because it 
is possible to go to large enough volume so that it is still numerically measurable and at 
the same time dominates over the subleading corrections by several orders of magnitude. 



6 Classical limit of scaling functions 



Motivated by the success of Mussardo et al. [11] who showed that in the case of periodic 
boundary conditions the (semi)classical limit of Liischer's finite size correction to the sine- 
Gordon soliton mass coincides with the infrared limit (/ = mL oo) of the soliton's 
classical energy in finite volume, we show next that the classical limit (s) of the previous 
Liischer type corrections to the ground state energy in the boundary sine-Gordon theory 
on the strip coincide with the asymptotic corrections to the classical energies of the various 
static (ground state) solutions of the model. 

In refs. ^ i4j collectively the following static solutions are described in details in terms of 
Jacobi elliptic functions: 

• The ground state (+) and the solitonic excited solution (-) with the $(0) = ^{L) = 
Dirichlet-Dirichlet (DD) boundary conditions where in the (+) [(-)] case the sine-Gordon 
field $(x) decreases [increases] monotonically from $^ at x = to [27r//3] at x = L; 

• the ground state (A) and the first excited state (B) type solutions in the zero topological 
charge sector with the 

DD boundary conditions, where, for the case (A) [(B)], ^(x) decreases [increases] from $^ 
to some turning point e [^ — e] then increases [decreases] from these values to 

• the ground state (-1) of the n = —1 topological charge sector with the 

TT /3$(L) 

< <^o < 2 ' — 2 — = '^■^ ~ ^ 
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DD boundary conditions, when decreases monotonically from positive to negative 

values, 

• and finally the ground state with mixed boundary conditions (DN): Dirichlet at x = 
and a one parametric (Ml) perturbed Neumann type b.c. at a; = L 

= ¥>o, f $'(^)U=L = -m^-^sin(^^), = 

For all these static solutions the dimensionless width of the strip I — mL and the 
energies of the solutions can be expressed in terms of the integrals 

b b 

l(a,b,C)= f , 7(a,6,C)= / dv\/sai^v^C. 

J Vsin' v + C J 

a a 

Indeed in the = Dirichlet-Dirichlet case for the (+) solution one obtains 



4?77 

I = 1(0, ^0, C), E+i^o, I) = -^m, Vo, C) - Cl/2], 



while for the (-) solution 



AfY) 

I = TT, C), E_{^o, I) = -^[1(^0, TT, C) - Cl/2] 



results. In the case of general DD boundary conditions, in the n = sector, for the type 
(A) (ground state) solution 

^ = XI ~ ^^"^ Ea{<Po: VL: I) = ^[ ^ -?"(e, Vi, - sin^ e) + ^ sin^ e] 

i=0, L ^ 1=0, L 

is obtained, while for the type (B) (excited state) solution 

I = XI ^ - - sin^ Eb{(Po: <PL: I) = ^[ X ^^^'^ TT - e, - sin^ e) + ^ sin^ e] 

i=0, L ^ i=0, L 

is found. In the n — —1 sector of the general DD boundaries for the only ground state we 
obtain 

/ = 1{Q, ipo, C)+l{Q, n-ipL, C), E^-'^ ((^o, ^l, I) = |^[/(0, <^o, C)+/(0, tt-^l, C)-Cl/2]. 
Finally for the ground state of the DN problem 

I = Vo, C), £;(^^)(¥^o, A, I) = ^[-C^- + n^L, <^o, C) - A~' cos(^l)] 
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is found with C = —(1 — A^"^) siri^ ipL- These expressions give the / dependence of the 
energies in a parametric form: in all cases the integration constant (s) C or e and (pi are 
- at least in principle - determined as functions of / (and the boundary parameters) from 
the first equations and using them in the second equations gives the / dependence of the 
energy expressions. 

The asymptotic form of the energies is obtained in the following way: first in each case 
the integrals defining I are expressed as an appropriate linear combination of complete 
and incomplete elliptic integrals. The / — > oo limit corresponds to the limit when the 
parameter of these elliptic integrals goes to unity - sometimes this is also correlated with 
the amplitude of the incomplete one going to 7r/2. Using the relevant asymptotic expansions 
of the elliptic integrals we determine the asymptotic C{1) (e(/), ipl)- In the second step the 
integrals appearing in the expressions for the energy are determined by realizing that their 
derivatives with respect to C (e, (p^) can be related to /(a, b, C); viewing this relation as a 
differential equation we integrate it (the constant is fixed by the known value of I {a, b,Oj). 
In this way the following expressions are obtained at / — cxd: 

E+ (y^o, /) ~ ^ (l - cos ipo + 4: tan^ ^ e"'^ 

E- (v^o, /) ~ ^ (l + cos (^0 + 4 cot' ^ e-^' 

Am ( LfQ 
Ea Wq, V^l, /) ~ (^2 - cos - cos v^l - 8 tan — tan — e" 

Eb (v^o, V^L, ~ ^ (2 + cos v?o + cos - 8 cot ^ cot 

,^^ 4m / iff) H — ifl _ 

E^ ' [ifo, V^L, /) ~ I 2 - cos cpo + cos cpi + S tan — tan e 

E(^^)(^o,^, - I? (1 - cos^o - - 4tan' YYT^^'~'') ' ^^'^^ 

The constant (/-independent) terms of these energy expressions are the classical energies of 
the static solitons/antisolitons forming the asymptotic ground state (for more details see 
[nmi). Note that some of the asymptotic corrections decrease with e~', while the others 
with e~''. The fact that the asymptotic correction is negative in the n = and positive in 
the n = —1 sectors can be understood by realizing that in the former sector the asymptotic 
solution is a superposition of a soliton and an antisoliton (whose interaction is attractive) 
while in the latter it is a superposition of two solitons (whose interaction is repulsive). 
The sign and magnitude of these asymptotic corrections were also verified numerically by 
comparing to the exact solution of the classical field equations. 

We now proceed to show that these asymptotic corrections are indeed the classical 
limits of the appropriate generalized Liischer formulae. 

It is very important to realize that the classical limit is not identical to taking 

/? ^ 0, in which almost all interesting quantities diverge (this limit cannot even be taken in 
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the classical formulae for the ground state energies). The general form of the semiclassical 
expansion for the ground state energy can be obtained from the Euclidean path integral 
as the limit 

If f 1 /'^/2 

E{L) = Jim -—log / D^exp <( -- / dr Lsi^) 



T-.00 hT 

where the Euclidean Lagrangian is 



h 



T/2 



+ Ml 



l-cos^($(L,r)-<l>i) 



+ Mo 



P 



l-cos|(<l>(0,r)-<l>o) 



Scaling out the coupling constant and introducing dimensionless variables 

$ = /?$ , X = mx , f = rriT 
standard semiclassical expansion around the classical solution yields 

Mo/32 Mi/52 



E{L) 



m 



F ml, 



m 



m 



(6.2) 



where F is a dimensionless function of its arguments. The first term is just the classical 
energy, the large volume asymptotics of which is given in (16.11) and is proportional to /3~2. 
All the quantum corrections are regular as /3 — 0, and expansion in powers ofh corresponds 
to expansion in P'^. As a result, the classical limit corresponds to isolating the leading 
term in the expansion around /3 = 0. Thus it is only meaningful to take the large volume 
asymptotic first and perform the classical limit by considering only the leading terms for 
small p. 

The volume dependence of the finite size correction is determined by the type of the 
particle giving the leading contribution: for terms of the form it is the first breather, 
while in the case of e"^' the second (note that m2 = 2mi = 2m in the classical limit). 

The UV-IR relations ()4.;^jl imply that the Dirichlet boundary condition is obtained as 
the ^ 00 limit of the general one, while the perturbed Neumann BC as the = 
{A~^ < 1) or 7] = {A~^ > 1) case. Using the classical UV-IR relation ()4.6p one finds that 
in the Dirichlet case 



Vcl,L = 



while for the perturbed Neumann BC 

COST]d=A~^ {A''^<1] 



coshT^ci = A~^ {A ^ > 1) 



results. 
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Using the above identifications, the classical limit of the breather one-particle term 
()4.9|) becomes 

32"^ tan ^ tan ^e-^ 

LcLli IcLli 1: 



/32 2 2 

which indeed agrees with Ea- 

To interpret Eb that belongs to the first excited boundary state let us recall that the 
refiection amplitudes on it can be obtained by the substitution 

7] — > f] = 71 {\ + 1) — T] (6.3) 

from the ground state ones [13 HH]. Therefore making this substitution in ()4.9j) gives the 
Liischer type correction to the first excited state energies. Since in the classical limit this 
substitution simplifies to rjd ^ — Vc\, it is clear that Eb does coincide with the classical 
limit of ()4.9|1 . Similar reasoning can be used in the case of E^-^K 

For the corrections of order e"^', however, things are more complicated. Let us take 
the example of the other cases can be treated in the same way. Taking the classical 
limit of the second breather one-particle term ()4.12ll gives 

32m 2 "Po -21 la a\ 

-^tan ye (6.4) 

which is off by a factor of 2. Fortunately this is not the end of the story. As we already 
noted above, in the classical limit the masses of the first and second breathers are given by 

mi = m , m2 = 2m 

and therefore the first breather two-particle term is exactly of same order as the second 
breather one-particle term. Let us examine the general formula for the two-particle con- 
tribution for a particle of mass m ((T 



rycosh0A^„(^^)K^(^^)e-2'"^^°^^^ = -me-2-^ T ^ /(^)e-2-^^'/2 
where we introduced 

f{e) = COshOK^ (6) KfS (6) Q-'^rnLicosh 9-1-9^2) 

Since f{9) is an even function, it expands in even powers of 9 around 

f{0) = h + he' + ... 

and so the two-particle contribution can be written as (for mL — > 0) 



aV^^ y ^ (2mL)" 
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At first sight this does not appear to be enough to rescue us, since the two-particle correc- 
tion is further suppressed by a factor of However, it is an interesting and nontrivial 
issue whether the asymptotic expansion (j6.5jl can be trusted in the semiclassical domain. 
The relevant integral in our case is 

/CO JQ poo jn 

— cosh^log(l + x(^)e-2"^^™^^'') ~ -mi / — cosh W e"'™^'^™''^' 
-oo 4vr 4:71 

(6.6) 



where 



and for the case of the boundary parameters are ?7o = v^o (-^ + 1), ^71, = and = 
{}l = oo. 

Using the known refiection factors, the coefficients fi can be explicitely calculated to 
have the following behaviour for small (3 



P4{n+l) 



which means that in the classical limit the series expansion ()6.5jl is invalid. 

The reason why the coefficients /„ diverge and the integral (I6.6jl must be treated very 
carefully is that the value of x(^) and all of its derivatives at 6* = go to infinity in the 
classical limit. The dangerous contribution comes from the particular block 

3 1 
2 ^ 2A 

in the refiection factor ()4.7|) and can be isolated from x as follows 

, , cosh 9 + cos , , 

X (0) = —r-. f Xo (0) 

cosh — cos ^ 

where Xo is a function with a smooth classical limit at the origin 
while the function 

, ^ cosh 6 + cos 

g{e) - 2^ 



cosh 9 — cos ^ 



has the Taylor coefficients 
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for large A. As a result, the series ()6.5j) takes the following form for small (3 



(l + f;(2n-l)!!( 

\ n=l ^ 




which is a power series with coefficients that grow factorially. It is also obvious that it does 
not have the right (3 dependence displayed in eqn. (|6.2|l . and the higher terms show more 
and more singular behaviour at /3 = 0. 

In order to get the correct asymptotics, we now evaluate the leading behaviour of the 
integral (|6.6p in the semiclassical regime carefully. First we substitute x = sinh6' for the 
integration variable. For large L, the integrand is concentrated around the origin and 
therefore we can evaluate the dominant contribution by keeping only the leading terms 
around a; = 0: 

1 - - A , ^ 1 + cos t + ^V2 . 



where a subtraction of —1 in the coefficient of Xo(0)e~^™^^ was applied to make the integral 
convergent at the plus and minus infinities (this only affects subleading terms - of order 
- as the singular behaviour around a; ~ is unchanged) . Now the integral can be evaluated 
using 



-nil 



TT TT TT 

sin^ — + cos ^Xo(0)e-2™i^ - sin — 



-^Xo(0)^e--.- 



4A 



using that L is very large. 

Taking the classical limit of small P (i.e. large A) and keeping only the leading term in 
P we obtain 

''^tan^^e--- 



/32 2 

which must be added to the second breather one-particle contribution ()6.4p . thereby giving 
perfect agreement with the classical formula for The same argument works for E^^^\ 
while -E_ can be obtained from by the substitution (jHSI)- 



7 Conclusions 

In this paper we succeeded in deriving the boundary analogue of Liischer's analytic for- 
mulae for finite size corrections (which were originally obtained for periodic boundary 
conditions). The consistency of these boundary Liischer formulae was checked against 
boundary TBA, NLIE and classical field theoretic results. The matching of all these dif- 
ferent approaches provided ample evidence for the correctness of the formulae, and also a 
very strong argument for the universality of the relation 
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between the one-particle coupling ga in the reflection factors and the one-particle amplitude 
Qa in the corresponding boundary state, which was flrst proposed for the boundary Lee- 
Yang model in [8] on the basis of matching the form factor expansion of one-point function 
against truncated conformal space data. Our numerical data based on the Dirichlet NLIE 
and the matching with the classical flnite size corrections show that this relation is valid in 
boundary sine-Gordon theory for very general values of parameters. The argument based 
on the boundary TBA shows that it is also valid for any theory with a diagonal bulk 
and boundary scattering theory. Based on these results we conjecture that this relation is 
indeed universal. 

Using the relation (IT.ljl . the boundary Liischer formulae proposed in this paper make 
it possible to compute the leading finite size corrections based only on data from the 
scattering amplitudes, even in non-integrable theories. It is interesting to note a deep 
analogy between the boundary Liischer formulae ()2.3|2.1flj) and the particle mass correction 
(11.111 for the periodic case. In both cases the generic correction can be written in terms 
of an integral over the relevant scattering amplitude, and the two integrals share many 
other similarities in structure. Besides that, the possible pole terms are also very similar. 
Therefore we expect that a quantum field theoretic derivation of our results along the lines 
of must exist. Such a derivation would be very interesting because it would permit 
generalization to any space-time dimensions (similarly to the case of periodic boundary 
condition while the cluster expansion technique is restricted at the moment to 1 -|- 1 
dimensional field theories, albeit not necessarily integrable ones. 

While the exponential dependence of the finite size corrections on the volume and the 
particle masses is obvious from general principles of field theory and TBA/NLIE, there are 
some highly nontrivial aspects of the results obtained here which we would like to stress. 
First of all, the relation (IT.ljl which was originally a (numerical) observation made in a 
particular model (Lee- Yang) in |S|, is shown to be universal by several pieces of evidence 
(analytic and numeric). Second, we now have a universal analytic form of the exponential 
finite size corrections expressed in terms of the infinite volume scattering theory, and also 
a regularized form for the two-particle contribution (eqn. (I3.7jl ) when the naive expression 
(I2.3jl is divergent. Third, we have done extensive consistency checks of the results, and in 
particular we have shown that they agree with classical field theory, but the classical limit 
of finite size corrections must be treated very carefully. 

In more practical terms, the boundary Liischer formulae are very useful for extracting 
boundary scattering data from finite size effects obtained by any numerical or analytic 
method like lattice field theory, truncated conformal space or Bethe Ansatz based ap- 
proaches (TBA, NLIE). Alternatively, when the scattering theory is explicitely known 
they can be applied to test the validity of numerical methods or conjectured TBA or NLIE 
equations. Such applications are presently being actively pursued by the authors. 
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00 


0L 


a 


^measured 


Tcxact 


A = 7/2 






exact: 0.867767 








0.05 


0.05 


0.867774 


-0.024529 


-0.024526 




0.03 


0.07 


0.867812 


-0.020640 


-0.020625 




0.07 


0.07 


0.867786 


-0.048403 


-0.048389 


A = 5/2 






exact: 1.175570 








0.1 


0.2 


1.176765 


-0.197730 


-0.193904 




0.05 


0.15 


1.175572 


-0.069682 


-0.069681 




-0.05 


0.15 


1.175571 


0.069681 


0.069681 



Table A.l: Comparing the first breather one-particle term to NLIE. For each value of A, 
we quote the exact value of a = rrii/M = 2sin^, which is compared to the exponent 
extracted from the NLIE. 7measurcd is the coefficient extracted from NLIE, while 7exact is 
the one predicted by (jA]lJ. 
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A Numerical comparison to the NLIE 

We consider the first breather one-particle contribution first. We normalize the energies 
and distances in units of the soliton mass M: 

EasiL) TT 1 + cos - sin Vo VL 9«ir, / /. x 

' = 2 sin— ^ ^ tan tan -^6-2"""^' + ... A.l 

M 2Al-cos§ + sin|^ 2A 2A ^ ' 

I = ML 
and the boundary values of the field as 

$. = ^0. , = -(-l)^/^0.7r (1 + A) 

Numerical data are obtained by fitting the NLIE results to a curve of the form 

76"°' 

Illustrative examples are given in table I A.l] the agreement was tested for many other values 
of the parameters involved, with results similar to those in the table. 

The second breather one-particle term proves more tricky. To eliminate the otherwise 
dominant first breather one-particle term we choose 0l = 0. However, the two-particle term 
for the first breather gives a non-negligible contribution for all values of the parameters that 
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10 


12 


14 


16 


A = 7/2, 00 = 0.1 


2 Bi 


-4.20321e-10 


-1.17697e-ll 


-3.35392e-13 


-9.67996e-15 




B2 


-1.87457e-09 


-8.21717e-ll 


-3.60199e-12 


-1.57893e-13 




E 


-2.29490e-09 


-9.39414e-ll 


-3.93738e-12 


-1.67572e-13 




NLIE 


-2.29710e-09 


-9.39742e-ll 


-3.93787e-12 


-1.67577e-13 


A = 9/2,00 = 0.013 


2 Bi 


-6.33781e-09 


-3.37019e-10 


-1.84416e-ll 


-1.03073e-12 




B2 


-1.31031e-07 


-1.00170e-08 


-7.65776e-10 


-5.85417e-ll 




E 


-1.37369e-07 


-1.03540e-08 


-7.84218e-10 


-5.95724e-ll 




NLIE 


-1.37425e-07 


-1.03550e-08 


-7.84232e-10 


-5.95718e-ll 


A = 9/2,00 = 0.007 


2 Bi 


-1.53201e-08 


-8.84074e-10 


-5.20261e-ll 


-3.10504e-12 




B2 


-5.51243e-08 


-4.21411e-09 


-3.22159e-10 


-2.46282e-ll 




E 


-7.04444e-08 


-5.09819e-09 


-3.74185e-10 


-2.77333e-ll 




NLIE 


-7.04671e-08 


-5.09861e-09 


-3.74193e-10 


-2.77334e-ll 



Table A. 2: NLIE versus second breather one-particle term combined with the first breather 
two-particle term 



are accessible by iterating the NLIE. Therefore we can only compare the NLIE results to 
the sum of the two terms. The first breather two-particle term can be safely calculated from 
(I2.3jl since the integrand is regular. Table |A]2l summarizes some of the results obtained this 
way. The agreement between the NLIE and the boundary Liischer formula becomes worse 
with decreasing volume, which is easy to understand since higher order terms of the cluster 
expansion start playing an increasingly important role. Once again, we have performed a 
rather extensive numerical check, of which the table contains only a small sample. 
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